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ABSTRACT 

The  mean  and  variance  of  the  product  of  the  narrowband  responses  of  two 
arrays  steered  to  the  same  look  direction  is  evaluated  and  shown  to  depend  on 
the  complex  coherence  between  the  two  array  outputs.  For  a  cross-line  array 
of  perpendicul ar  equi-spaced  elements,  the  stability  is  much  poorer  than  for  a 
sum  array  processor,  due  to  the  largely  uncommon  volume  of  intersection  of  the 
two  cones  of  response  of  each  line.  When  each  array  is  extended  to  be  planar, 
the  degradation  in  stability  is  lessened,  tending  towards  the  sum  array 
performance  as  the  number  of  common  elements  increases. 
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INTRODUCTION 

A  planar  array  with  a  grid  structure  of  x  M2  elements  requires  a 
large  number  of  receivers  and  considerable  signal  processing  when  all  the 
elements  are  employed  and  actively  utilized.  In  an  effort  to  conserve  on  the 
number  of  elements  and  amount  of  signal  processing,  the  possibility  of  using  a 
sparse  array  seems  to  have  merit.  In  particular,  it  has  been  found  that  an 
equi-weighted  planar  array  has  the  same  auto  spectral  density  response  as  the 
cross  spectrum  of  a  pair  of  perpendicul ar  lines  of  double  the  length  and  with 
triangular  weighting  on  each  line.  However,  without  investigating  the 
variances  of  these  sum  and  cross-line  arrays,  respectively,  it  is  impossible 
to  decide  on  their  relative  merits. 

Here  we  will  consider  two  arbitrary  planar  arrays  (which  may  have  some 
common  elements  and  may  even  be  linear  arrays),  both  of  which  are  steered  to 
the  same  look  direction  and  employ  weightings  for  sidelobe  control.  The 
sample  cross-spectral  density  of  the  two  array  outputs  is  the  output  variable 
of  interest.  Thus  the  combination  of  arrays  is  resolving  in  both  spatial 
angles  (wavenumber)  as  well  as  in  temporal  frequency.  We  will  evaluate  the 
mean  and  variance  of  the  cross-spectral  density  estimate  at  the  system  output 
in  terms  of  the  statistical  properties  of  the  impingent  noise  field  and  the 
array  parameters,  such  as  look  direction  and  weighting. 

As  a  special  case,  by  choosing  the  two  planar  arrays  identical  in  element 
usage  and  weighting,  we  will  reduce  to  the  sum  array  since  the  cross-spectrum 
then  becomes  the  auto-spectrum.  Thus  we  can  compare  the  performances  of 
product  arrays  and  sum  arrays  in  terms  of  the  mean  and  variance  of  their 
responses. 
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CROSS-SPECTRAL  DENSITY  ESTIMATE 

Let  x(t)  and  y(t)  be  any  two  array  outputs  which  are  stationary  in  time, 
zero-mean,  and  have  auto-spectra  and  cross-spectrum  Gx(f),  Gy(f),  GXy(f), 
respectively,  f  is  temporal  frequency  in  Hz.  Sections  of  each  waveform  are 
gated  out  by  multiplying  by  temporal  weightings  and  then  subjected  to  Fourier 
analysis  according  to* 

Xm(f )  =  Jdt  exp(-i2irft)  wm(t)  x(t)  for  1  <_  m  <  N, 

Yn(f)  =J'dt  exp(-i2irft)  wn(t)  y(t)  for  1  <_  n  <_  N.  (1) 

These  temporal  weightings  are  generally  taken  as  delayed  versions  of  a  basic 
weighting  w(t)  according  to 


wn(t)  =  w(t-nS) ,  (2) 

where  S  is  a  shift  or  time  delay;  however,  we  keep  the  more  general  case  in 
(1)  for  the  time  being. 

The  cross-spectral  density  estimate  at  frequency  f  is  obtained  by 
multiplying  outputs  (1)  and  averaging  in  time  according  to 

N 

v-  5  Vf>  Vn(f)'  <3) 

n=l 


This  is  the  product  array  output;  it  has  resolution  capability  in  the  spatial 
angles  by  virtue  of  each  array  output  being  steered  to  the  same  desired  look 
direction,  and  it  has  frequency  resolution  governed  by  the  lengths  of  the 
weightings  in  (1)  or  (2). 

The  mean  of  a  general  product  term  of  components  of  (1)  is 


*An  integral  without  limits  is  over  the  entire  range  of  its  non-zero  integrand. 
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Xn/f^Vn^  “  JTdt  du  exp[-12.f (t-u)]  wm(t)  wn(u)  Rxy(t-u)  = 

-  |'<«-exp(-12,fr)  Rxy(T)  im(T)  -  8xy(f)«fmn(f)  = 

*  [du  Gxy(u)imn(f-u)*  (4) 

where 

RxyCC)  =  x(t)  y(t-T)  (5) 

is  the  cross-correlation  of  array  outputs, 

fmn<T>  -  Jdt  wn(t-D  (6) 

is  the  aperiodic  correlation  of  wm  and  wn,  ®  denotes  convolution,  and 

"  j<H-exp(-i2*fT)  *mM  -  "ra(f)  «*(f)  (7) 

is  the  product  of  windows  of  the  individual  temporal  weightings. 

If  cross-spectrum  G  does  not  vary  significantly  in  the  width  of 

X  J 

window  $mn,  (4)  yields  approximation 


X  (f)Y  (f)  a  G  (f) 
nr  '  nv  xyv  ' 


r  _ 

du  ©  (f-u)  =  G  (f)  0  (0). 

*mnv  '  xyv  ’  nnnv  ’ 


(8) 


Now  we  apply  these  results  to  find  the  mean  of  the  system  output  cross-spectral 
density  estimate  in  (3): 


-  5  -  Sxy<f) 


N 

5 

n=l 


0  (0)  =  G  (f) 

^nnv  '  xyv  ' 


N  r 

z* 

n=l 


dt  wn(t). 


(9) 


which  is  proportional  to  the  true  cross-spectral  density  G  between  the  two 

xy 

array  outputs. 
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In  order  to  evaluate  the  variance  of  system  output  v  in  (3),  we  need 
magnitude-square  value  (suppressing  f) 


N 


m,n=l 


\n  Ym  Xn  V 


(10) 


We  assume  that  the  filtered  outputs  Xm  and  Yn  in  (1)  are  complex  Gaussian 
(as,  for  example,  if  x(t)  and  y(t)  were  Gaussian).  We  then  observe  that 


XmYn  =  jTdt  du  exPC-i2irf (t+u)]  wm(t)  wn(u)  Rxy(t-u)  = 

=  Jdv  Gxy(v)  Wm(f-v)Wn(f+v).  (11) 

Now  for  f  removed  from  zero  by  at  least  the  reciprocal  of  the  segment  length, 
the  two  windows  in  (11)  do  not  overlap,  and  we  get 


Xm  Yn  a  0  for  f  *  0  . 


(12) 


Then  using  the  factoring  property  of  zero-mean  Gaussian  random  variables, 
(10)  becomes 


m 


N  r _ „ 

j>  fx  Y* 
^  I  m  m 

,n=l  L 


X  Y*  X*  Y  +  Xm  Xn  Ym  Y 
m  m  n  n  m  n  m 


•k  k 

X  Xn  Y  Y  , 
m  n  m  n  ’ 

— 


(13) 


where  we  used  (12).  The  variance  of  random  variable  v  is 

N 


°v  =  M7  =  M  ^  -  |  v  |2  = 


Xm  Xn  Ym  Yn 


m,n=l 


>  i 

m,n=l 


[Gx(f)«imn(f)][Gy(f)®$nln(f)] 


=  s  (f)  G  (f)  'y  i  2(o) 
x  y  X  mn 

m,n=l 


(14) 
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where  we  used  (9),  (4),  and  (8). 


We  now  specialize  the  general  weightings  in  (1)  to  the  particular  case  in 
(2) ,  obtaining  from  (6) 


where 


dmn(0)  =  fdt  w(t-mS)  w(t-nS)  =  dw((n-m)S), 


tfw(T)  =  J'dt  w(t)  w ( t— t) 


(15) 


(16) 


is  the  aperiodic  correlation  of  basic  temporal  weighting  w.  Then  variance 
2 

av  in  (14)  becomes 


N  N-l 

al  =  Gx(f)  Gy(f)  ^^((n-m)S)  =  Gx(f)  Gy(f)  (N-jpj)  «^(pS)  .  (17) 

m,n=l  p=l-N 

We  are  now  in  position  to  formulate  a  quality  ratio  for  the  output  of  the 
product  array.  Namely  we  define  the  complex  (voltage)  quality  ratio 


Y  (f ) 
rxy\ 


(i 


JSL) 


<*w(pS) 


__-l/2 
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where  we  define  the  complex  coherence  between  the  array  outputs  as 


(18) 


Yxy 


(f) 


Wf ) 

[Gx(f )  Gy(f)]1/Z  ’ 


(19) 


and  have  employed  (9),  (17),  and  (15).  The  quantity  in  (18)  is  desired  large; 
it  has  leading  factor  yif,  which  however  is  partially  compensated  by  the  last 
factor  of  (18)  if  shift  S  is  less  than  the  segment  length  of  temporal  weighting 
w.  A  detailed  investigation  of  the  temporal  processing  factors  in  (18)  is 
given  in  [1];  for  present  purposes,  the  factor  is  nearly  maximized  if  shift  S 
is  taken  about  50% of  the  segment  length. 
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However,  the  most  important  quantity  in  quality  ratio  (18)  is  the  complex 
coherence  at  frequency  f,  y  (f).  It  is  always  bounded  in  magnitude  by  1, 
and  can  be  significantly  less  than  1  if  the  two  array  outputs  x(t)  and  y(t) 
are  incoherent  at  frequency  f  of  interest. 

On  the  other  hand,  for  identical  arrays  and  array  weightings,  we  have 
y(t)  =  x(t),  and  the  quality  ratio  is  again  given  by  (18),  where  YXy(f)  is 
replaced  by 


Yxx(f)  -  1. 


(20) 


That  is,  the  sum  array  and  product  array  differ  in  their  complex  quality 
ratios  simply  by  the  factor  vXy(f),  which  is  the  complex  coherence  of  the 
two  arrays  in  the  product  formulation.  Thus  the  relative  performance  of  a 
product  array  can  be  investigated  by  determining  the  coherence  of  its  component 
array  outputs. 

With  this  information,  we  can  now  give  a  qualitative  measure  of 
performance  of  the  cross-line  array.  Consider  two  line  arrays  lying  along  the 
horizontal  and  vertical  axes,  respectively.  Suppose  both  lines  are  steered  to 
the  same  look  direction  in  three.-dimensional  space.  Since  a  line  array  must 
inherently  have  maximum  response  everywhere  in  a  cone  of  symmetry  centered  on 
the  line,  the  two  cones  will  intersect  at  the  desired  look  direction,  but  will 
have  largely  non-overlapping  cones  at  other  angles.  Thus  only  a  small  fraction 
of  the  output  of  each  array  is  in  common;  in  fact,  most  of  each  array  output 
comes  from  uncommon  arrival  angles. 

Thus  the  two  line  array  outputs  will  be  largely  independent  of  each 
other,  meaning  low  coherence.  Furthermore,  the  longer  the  line  arrays,  the 
finer  becomes  the  angular  resolution,  and  the  common  intersection  volume  of 
the  cones  decreases.  Thus  the  performance  of  the  product  array  relative  to 
the  sum  array  becomes  poorer  as  the  line  arrays  become  longer.  In  the  next 
section,  these  conclusions  will  be  verified  by  a  detailed  quantitative 
analysis  of  the  coherence  between  two  general  array  outputs. 
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CROSS-SPECTRUM  OF  INDIVIDUAL  ARRAY  OUTPUTS 

Let  the  pressure  field  at  time  t  and  general  location  x,y  in  a  planar 
array  be  denoted  by  p(t,x,y).  Let  the  field  be  stationary  and  homogeneous, 
with  temporal-spatial  correlation 


p(tpX^y^)  pTt^Tx^Ty^T  —  Rp(t^  1 2 * x^— x 2 ^1” y 2. )  *  (21) 

The  frequency-wavenumber  spectrum  corresponding  to  Rp  is  <j>p(f,u,v),  where 

Rp(T,u,v)  =  JJ/df  du  dv  exp( i2irfr+iuM+ivv)  <^(f,u,v)  .  (22) 

We  also  define  a  partial  transform  of  (22)  as 

J0fp(f.u,v)  =  JJdu  dv  exp(  iuM+ivv)  5p(f  ,u,v)  .  (23) 

This  mixed  function  of  temporal  frequency  and  spatial  separations  will  be  of 
prime  importance  later,  especially  if  it  can  be  evaluated  in  closed  form. 

The  grid  structure  of  the  planar  array  is  such  that  the  elements  are 
equi-spaced,  being  located  at  positions  md,nd  in  the  x,y  plane,  for  m,n 
integer.  If  a  particular  element  is  absent  or  is  not  used  in  an  array  output, 
the  weighting  of  that  element  output  is  simply  set  equal  to  zero.  For  polar 
angle  <t>  as  measured  from  the  z-axis,  and  azimuthal  angle  e  as  measured  from 
the  x-axis,  the  time  delay  employed  at  location  md,nd,  in  order  to 
steer  in  desired  look  direction  ,is  [2,  Appendix  A] 

rmn  “  “  £(am+6n)  .  (24) 

where  d  is  the  element  spacing,  c  is  the  speed  of  propagation,  and 

a  =  si n<t>ji  coseg  ,  b  =  sin^  sine^  .  (25) 

Broadside  to  the  planar  array  corresponds  to  ^  =0. 
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Array  output  x(t)  is  synthesized  by  choosing  a  particular  subset  of 

elements  in  the  planar  grid  structure,  weighting  their  outputs,  and  time-delay 

* 

steering  to  the  desired  look  direction  ^  .  Thus 

x(t)  =^>wx(k’^  p »  kd>  Jld)  .  (26) 

W 

By  simply  setting  some  weights  to  zero,  a  line  array  or  a  cross-line,  or  any 
desired  array  configuration,  can  be  realized.  Similarly,  a  second  array 
output,  which  may  employ  some  or  all  of  the  same  elements,  is 

y(t)  =^wy(m’n)  md>  nd)  ‘  (27) 

mn 

The  cross-correl ation  of  array  outputs  (26)  and  (27)  is,  upon  use  of  (21), 

Rxy(T)  =  x(t)y(t-T)  =  wx(k,^)  wy(m,n)  Rp(r+T  n-T^  ,  (k-m)d,  U»-n)d) .  (28) 

Wmn 

Using  (23)  and  (24),  the  cross-spectrum  of  x  and  y  is 

Gxy(f)  =  wx(kj)  wy(m,n)  exp[i2irf(rnn-ZjJJl)]^p(f,(k-m)d,CiL-n)d)  = 
k^mn 

=  ]!>  wx(kj)  w  (m,n)  exp[i2irf^-(ak+6^  -om-Bn)]  Jrp(f , (k-m)d, (/-n)d)  = 
iy£mn 


=  ^^xy(q,r)  exp|l2irf|(oq+BriJ  ^p(f,qd,rd)  , 


(29) 


qr 


where  we  let  q=k-m,  r=  -n,  and  defined 


xy' 


tfyw(q.r)  wx(kj)  w  (k-q,^-r) 


(30) 


k* 


*A  summation  without  limits  is  over  the  entire  range  of  its  non-zero  summand. 
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as  the  two-dimensional  cross-correlation  of  the  weight  structure  employed  to 
yield  array  outputs  x(t)  and  y(t). 

As  the  first  special  case  of  the  above,  a  pair  of  perpendicular  cross- 
lines  can  be  realized  by  setting 

wx(k,^)  =  0  except  for  £=  0, 

wy(m,n)  =  0  except  for  m  =  0;  (31) 

in  this  case,  (30)  yields 

tfxy(q,r)  =  wx(q,0)  wy(0,-r),  (32) 

and  (29)  becomes  the  cross-spectrum  of  the  cross-line  array, 

Gxy^(f)  =5WX((1’0)  wy(°’~r)  exp[i2irf|-(aq+8r)]^)(f,qd,rd).  (33) 

qr 

A  second  special  case  is  obtained  by  considering  the  sum  array,  which 
corresponds  to  setting 

wx(k,i)  =  wy(k,Jl)  =  w(kj), 

thereby  getting  from  (30) 

tfxy(q,r)  =^w(k,^)  w(k-qU-r)  s  ds(q,r)  , 
k£ 

and  from  (29),  the  auto-spectrum  of  the  sum  array, 

Gxy^(f)  =^^s(q’r)  exp[12irf^(aq+Br)]  Jrp(f,qd,rd)  .  (36) 

qr 

Comparison  of  special  cases  (33)  and  (36)  reveals  that  the  mean  responses 
of  the  cross-line  and  sum  arrays  can  be  made  equal  by  setting 


(34) 


(35) 
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wx(q,0)  wy(0,-r)  =  0s(q,r)  =2  w(k,/0  w(k-qj-r)  for  all  q,  r  .  (37) 

So,  for  example,  if  the  sum  weighting  w(k,^)  is  flat  over  a  rectangle, 

$s(q,r)  is  triangular  in  q  as  well  as  r,  over  a  rectangle  twice  the  size. 

Thus,  if  line  weightings  w v(q,0)  and  w..(0,-r)  are  each  triangular  over 

this  double-length,  (37)  is  satisfied,  and  Gx^(f)  =  G^(f).  This  conclusion 

holds  irrespective  of  the  pressure  field  statistics  J&p(f ,qd,rd) . 

Returning  now  to  the  general  case  of  (29)  and  (30),  the  auto-spectrum 
Gx(f)  of  array  output  x(t)  in  (26)  is  easily  obtained  by  replacing  w^  by 
wx  in  (30)  and  using  that  result  for  0xy  in  (29).  A  similar  procedure, 
but  now  replacing  wx  by  wy,  yields  the  auto-spectrum  Gy(f)  of  array 
output  y(t)  in  (27).  Combined  with  (29)  itself,  we  now  have  the  capability  of 
calculating  coherence  yx  (f)  as  given  earlier  by  (19). 

An  alternative  illuminating  form  for  cross-spectrum  G  (f )  is  obtained 
n 

by  substituting,  for  in  the  second  line  of  (29),  the  expression  (23)  and 
interchanging  summation  and  integration: 

GXy(f)  =  jjdy  dv  <£p(f ,y,v)* 

*W  ( 2Trf-^a+dy,  2ir4fl+dv)  W*(2ir4a+dy,  2ir4e+dv)  ,  (38) 

A  L.  y 


where 

Wx(a,b)  =  2  *X(M)  exp(ika+Ub), 

W  (a,b)  =  ^  wy(m»n)  exp( ima+i nb)  ,  (39) 

mn 

are  the  response  patterns  of  the  x  and  y  arrays.  Since  the  patterns  in  (39) 
peak  at  a  =  b  =  0,  the  integral  in  (38)  is  dominated  by  the  contribution  at 
spatial  frequencies 
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That  is,  the  cross-spectrum  in  (38)  is  influenced  mainly  by  frequency- 
wavenumber  spectrum  value  Jp(f,u,v)  at  the  values  given  by  (40).  x=c/f  is 
the  wavelength  at  the  frequency  f  of  interest. 
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EXAMPLES  OF  FREQUENCY-WAVENUMBER  SPECTRA 


Square  Support 


Suppose  that  at  some  frequency  f=fQ,  the  frequency-wavenumber  spectrum 
is  flat  over  a  square: 

il(fo)  ^2  for  W  <  kL’  W 

/  0  otherwise 

where  k^  is  the  common  limit  of  wavenumbers  in  p  and  v.  Then  from  (23), 

sin(k,  u)  sin(k.  v) 

<fo’u*v)  ■5l(fo)  "T[r - EjT^-  <«) 

This  closed-form  expression  is  separable  in  u  and  v  and  will  lead  to 
worthwhile  simplifications  when  employed  in  (29),  and  especially  cross-line 
result  (33). 


Circular  Support 


Suppose  instead  that 


fp(f0>w,v) 


(ll<fo>  “If  for  vi2  +  v2  <  k^ 

irk£ 


Then  (23)  yields 


0  otherwise 


ip(f0.u,v)  -  Si.(f0) 


k^  +  / 


(43) 


(44) 
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which  has  circular  symmetry  in  separation  space  u,v.  Although  in  closed  form, 
it  is  not  separable  in  u  and  v,  and  is  more  time-consuming  to  evaluate  than 
(42). 
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CROSS-LINE  ARRAY  WITH  SQUARE-SUPPORT  SPECTRUM 

In  this  section,  we  specialize  to  a  cross-line  array  as  depicted  in 
figure  1  and  couple  it  with  the  square-support  frequency-wavenumber  spectrum 


of  (41).  The  cross-line  array  has  one  element  at  the  origin  and  M  elements  in 
each  of  the  four  perpendicul ar  legs  protruding  from  the  origin,  for  a  total  of 
4M+1  elements.  All  elements  are  spaced  by  d  in  both  directions,  consistent 
with  ( 24)-( 27) . 

The  cross-spectrum  at  frequency  f  ,  between  the  two  perpendicul ar 
lines,  is  given  by  (33)  with  (42): 


Gxy^V  =?l<fo)  £wx(q>°)  wy(°»-r)  exp[ikQd(aq+sr)]  — 


sin(k^dq)  sin(k^dr) 


k^dr 


qr 


M  sin(k.dq)  M  sin(k.dr) 

=  4il(f„)2  «qwx(q,0)cos(.kodaq)  — jrgg - ^  «rwy(0,r)cos(kod8r)  -j-g; - ,(45) 

q=0  L  -  n  L 


r=0 


where 


k  =  —  =  2*f _/c  . 

1 

,1/2 

for 

q  =  Ol 

o 

o 

o 

'x 

for 

q  >  lj 

(46) 
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and  the  weight  structure  on  each  line  has  been  assumed  real  and  symmetric 
about  the  origin.  The  quantities  in  (45)  are  all  real  and  require  only  two 
single-summations  of  size  M  at  each  value  of  the  three  dimensionless 
parameters  kQda,  kQd8,  k ^d .  Here,  is  the  common  limit  on  allowed 
wavenumbers  in  square-support  spectrum  (41),  and  kQ  is  the  wavenumber 
corresponding  to  frequency  f  of  interest.  Quantities  a  and  8  are  given  in 
terms  of  the  look  direction  according  to  (25). 


To  find  the  auto-spectrum  of  line-array  output  x(t),  we  replace  w^  by 
wx  in  (30)  and  use  the  upper  line  of  (31): 


(47) 


where 

^xx(q)r  ]>wx(k,°)  wx(k-q,0)  (48) 

k 

is  the  auto-correl ation  of  the  x-array  weights.  Substitution  of  (47)  in  (29) 
yields 


GiC>(fo)  =  5.^xx(q)  exp(ikodaq)j2p(fo,qd,0)  = 

q 


2M  sin(k, dq) 

-  2  £l<fo>  j>  'q  c°s(kodaq)  v-fr—  (49) 

q=o 

where  we  used  (42)  and  (46).  In  a  similar  fashion,  the  auto-spectrum  of 
line-array  output  y(t)  is  given  by 

( c )  t-  ,  sin(k.dr) 

y  <V  =  2^(fo>  2  er  Vr)  cos(kod6r>  k.  dr  •  (50) 

r=0  L 


where 
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fyy(r)s]>wy(0,J?)  wy(0,^-r)  ,  (51) 

in  keeping  with  (31).  The  auto-spectral  results  in  (49)  and  (50)  each  require 
real  single-sums  of  size  2M.  No  double  summations  are  required  in  the  cross- 
spectral  result  of  (45)  or  in  the  auto-spectral  results  of  (49)  and  (50), 

although  in  the  latter  cases,  the  auto-correlations  l^x  and  must  be 

pre-computed.  The  complex  coherence  at  frequency  fQ,  y*  '  (fQ) ,  between 

individual  array  outputs  x ( t )  and  y(t)  of  the  cross-line,  is  given  by  ratio 
(19)  as  usual;  the  factor  ^(fQ),  as  well  as  the  absolute  scales  of  the 
weight  structures  (wx(k,0)}  ,  £wy(0,i)},  will  cancel  out  in  the  ratio.  The 
fundamental  parameters  are  kQda,  kQdB,  k^d. 

For  broadside  steering  of  the  cross-line  array,  we  have  ^  =  0  and  (25) 

yields  a  =  8  =  0.  Then  (45),  (49),  (50)  are  independent  of  kQ,  and  the 

coherence  depends  only  on  k^d. 
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CROSS-LINE  ARRAY  WITH  CIRCULAR-SUPPORT  SPECTRUM 

By  combining  (33)  with  the  circul ar-support  frequency-wavenumber  spectrum 
of  (44),  and  using  assumed  symmetry  of  the  line-array  weights,  we  obtain  cross 
spectrum 


.(c) 

axy 


M  M  2J; 

(f0)  =  4$i(fo)  eqwx(q,0)cos(kodotq)  erwy(0,r)cos ( kQd6r) 

q=0  r=0 


'1(kLd|/q2+r2  ) 


LdF 


2  2 
c+r 


(52) 


This  double  sum  can  no  longer  be  separated  into  the  product  of  two  single 
sums,  as  (45)  was,  due  to  the  coupling  caused  by  the  Bessel  function. 

However,  the  Bessel  function  need  only  be  computed  in  a  45°  sector  of  the  q,r 
plane  and  then  reflected  about  the  45°  line;  i.e.,  the  same  value  is  attained 
for  q,r  =  m,n  as  for  q,r  =  n,m. 


The  auto-spectrum  of  the  x-array  output  is  obtained  by  utilizing  (44)  in 
the  top  line  of  (49): 


fi‘C)(f0)  -  2fj(f0)  2.  eq  cos(k0daq) 

q=0  L 

where  ^x  is  again  given  by  (48).  In  a  similiar  fashion,  there  foil 


2M 


2J1(k1_dq) 


(53) 


ows 


Gyc)(f0)  =  2?1  (f0)  %  erfyy(r)  cos(k0der) 


2M 


2Jx(kLdr) 


y 


yy 


k,  dr 


(54) 


r=0 


The  coherence  is  now  available  from  ( 52)— ( 54 ) .  This  example  was  not  pursued 
numerical ly. 
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RESULTS 

A  program  for  the  calculation  of  the  coherence  of  a  cross-line  array,  via 
(45 )— (51 ) ,  is  given  in  appendix  A.  It  was  exercised  to  give  the  following 
results  in  figures  2-6.  Although  these  plots  look  virtually  identical,  closer 
inspection  of  the  numerical  values  (not  included)  reveals  that  there  are  small 
(insignificant)  differences,  even  for  the  widely  different  values  of  kQd, 

,  ©p  considered  here.  The  degradation  of  the  cross-line  array  relative  to 
the  sum  array  is  virtually  independent  of  the  particular  look  angle  ^  ^  . 

This  can  be  partially  explained  by  virtue  of  the  fact  that  at  broadside 
steering,  a  line  array  has  a  narrow  beam  but  covers  a  full  360°  angular  sweep, 
whereas  at  endfire,  the  beam  is  broad  but  exists  at  only  one  angle.  Thus  the 
total  angular  coverage  is  essentially  constant. 

The  overriding  impression  of  figures  2-6  is  that  the  degradation  of  the 
cross-line  array  is  significant,  relative  to  a  sum  array,  in  terms  of  the 
stability  of  the  cross-spectral  estimate.  This  is  particularly  so  as  the  size 
of  the  array  (4M+1  total  elements)  grows,  or  as  kLd  increases  above  .5. 

Further  examples  of  interest  may  be  obtained  from  the  program  in  appendix  A. 
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FUjvinr  3.  D^radaFion  -far  k0  d  =  ^  ?  ^  -  0 
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SUMMARY 


When  the  horizontal  line  array  in  figure  1  is  time-delay  steered  to  look 
in  some  desired  direction  ^ ,  it  must  also  respond  in  a  cone  of  symmetry 
centered  on  the  axis  of  the  line.  Similarily,  although  the  vertical  line  is 
steered  to  look  in  the  same  direction  ,  it  too  has  a  cone  of  equal 

response,  but  now  centered  on  the  vertical  axis.  These  two  cones  will 
intersect  at  and  thereby  lead  to  some  common  power  at  frequency  f  at 

their  respective  array  outputs,  prior  to  multiplication  and  averaging 
according  to  (3).  However,  both  arrays  also  respond  to  uncommon  (i.e. 
uncorrel ated)  power  contributions  from  other  directions,  each  within  its  own 
cone  of  response.  The  sharper  the  beams  of  each  line  array,  the  less  common 
power  will  be  intercepted,  leading  to  less  coherence  between  the  two  line 
outputs.  Thus  the  stability  of  the  cross-line  array,  relative  to  the  sum 
array,  degrades  as  the  size  of  each  line  array  increases. 


If  each  of  the  product  arrays  were  made  of  a  parallel  pair  of  lines 
(separated  by  some  integer  multiple  of  d),  the  responses  of  each  would  be 
rather  complicated.  However,  there  would  again  be  some  common  overlap  at 

,  but  a  great  deal  of  uncommon  response  at  other  sidelobe  regions,  still 
causing  a  decreased  coherence  and  unstable  estimates.  As  more  parallel  lines 
are  added  to  each  product  array,  the  performance  should  monotonical ly  approach 
that  of  the  sum  array,  being  actually  realized  when  all  of  the  available 
parallel  lines  are  employed,  since  the  outputs  of  the  horizontal  and  vertical 
arrays  are  then  identical. 
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APPENDIX  A.  PROGRAM  FOR  CROSS-LINE  ARRAY 

The  following  program  calculates  the  coherence  of  a  cross-line  array  by 
means  of  (45)— (51).  M  is  the  number  of  elements  in  each  of  the  four 
perpendicul ar  legs  protruding  from  the  origin.  The  weighting  in  line  60  is 
assumed  the  same  for  both  lines  and  is  triangular.  The  autocorrel ation  of  the 
weight  structure  is  computed  and  stored  in  line  130.  Inputs  of  kQd, 
are  required  in  lines  160-180  respectively.  When  these  latter  inputs  are 
desired  changed,  the  program  can  be  continued  at  line  160  instead  of  10, 
provided  that  the  array  weights  have  not  also  been  changed.  The  input  of 
k^d  occurs  in  line  320.  The  dB  output  in  line  520  is  according  to 

dB  =  10  log10jrxy(fo)|2  . 

since  y  is  proportional  to  the  voltage  quality  ratio. 


10 

D I M  W < 0 : 3 0 , - 3 0 : 3 0 > , Psi < 0 : 

3  0 , 0 :  6  0  > 

2  0 

D  I  ri  C  o  2-  a  <  8  :  6  0  >  ,  C  o  s  b  <  0  :  6  0 ) 

, S i nc <0 : 69) , 

,  Dbt 

!  6  :  3  0  i 

3  8 

FOR  M=0  TO  30 

40 

M  1  =  M+  1 

50 

POP  Ms  =  0  T 0  M 

60 

W  < M , M s >  =  W ( M ,  - M s )  =  1- M s  ■  -  M 1 

i  a  rt 

■■  ay 

w  e  i  g  h  t  i  n  g 

70 

NEXT  Ms 

S  0 

FOR  Q s  =  0  TO  M  +  M 

90 

S  =  0 

1  0  0 

FOP  Ms-Qs-M  TO  M 

1  10 

S  =  S  +  W  <  M  ,  M  s  >  *  W  <  M  ,  M  s  -  Q  s  > 

120 

NEXT  Ms 

130 

Psi  (M, Qs>=S  ! 

au  t  o  ~  c  o  r  r  t* 

1  at  ■ 

ion  ot'  anna 

140 

NEXT  Qs 

150 

NEXT  M 
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160  Kod  =  P  I  i  5  ! 

170  P o  1  ar  =  P  I  ••••' 4  ! 

1 3  0  fl  z  i  m  u  t  h  =  P  I  x  4  ! 

1  9  O  3  —  K  o  d  S  I  N  (.  P  o  1  ar  ) 

2 0 @  T  a=  3 * C O S <flzi  m u t h > 

210  Tb  =  S * S  I N  (.  fl z  i  rn u t  h  > 

2  2  0  F  0  R  Q  s  =  0  T  0  6  0 

2  3  0  C  o  s  a  (  Q  s  >  =  C  0  S  <  T  a  *  Q  s  ) 

2  4  0  C  o $ b < Q s ) =  C O  S ( Tb* Q  s ) 

250  NEXT  Qs 

260  PLOTTER  IS  "93720" 

270  L I M I T  3@ , 170,40, 235 
230  OUTPUT  735; " VS 2  " 

2  9  0  S  C  fl  L  E  0 , 3  0  ,  -  3  0 , 0 

3  0  0  GRID  5,5 
310  PEHUP 

320  FOR  K 1 d  = . 25  TO  3  STEP  .25 

3  3  0  S  i  n  c  l.  0  >  =  1 

3  4  y  F  U  R  0  3—1  T  u  6  y 

3  5  0  S  =  K  1  d  *  Q  s 

3  60  Si  n  c  C  Q  s  )  -  S  I N  <  S  >  /'  S 

370  NEXT  Q s 

3  3  0  F  0  R  M  =  0  T  0  3  0 

3  9  0  S  1  a=  S  1  b  =  ,  5  *  W  (  M  ,  0  > 

400  FOR  Q s  = 1  TO  M 

410  S  =  W <  M , Q  % ) * S i n  c < Q  s  > 

4  2  0  S  1  a  =  S  1  a  +  S  *  C  o  s  a  (  Q  s  > 

4  3  0  S  1  b  =  S  1  b  +  S  *  0  o  s  b  (  Q  s 
440  NEXT  Qs 

4 5 0  S  2  a=  S  2  b  =  .  5  *  P  s  i  <  M  ,  0  > 

4  6  @  F  Ci  R  Q  s  =  1  T  0  M  +  M 
4  7  0  S  =  P  s  i  <.  M  ,  Q  s  >  *  S  i  n  c  <  Q  s  ) 

4  8  0  S  2  a  =  S  2  a  +  S  *  C  o  s  a  (  Q  s  > 

4  9  0  S  2  b  =  S  2b  +  S  *  C  o  s  b < Q  s ) 

5  0  0  N  E  X  f  Q  s 

5  1 0  0  o  h  e  r  €•  n  c  e-  =  2  *  S  la*  S 1  b  S  Q  R  (  $ 

5  2  0  D  b  (  M  =  1  0  *  L  G  T  d  0  o  h  e  r  &  n  c  <s  A  2  . 

530  NEXT  M 

5  4  0  F  0  R  M  =  0  T  0  3  0 

550  PLOT  NjIib'id'O 

560  NEXT  M 

570  PEHUP 

580  NEXT  Kid 

590  END 


INPUT 

INPUT 

INPUT 
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